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Abstract 

Motivated by astrophysical observations of excess electronic production in the galaxy, 

■ new theoretical models of the dark matter sector have been proposed in which the Standard 

O ■ 

model couples to the dark matter sector through an attractive interaction. The coupling 



(N 



dt 



of the Standard model to "dark strings" , which are solutions of the low energy dark sector 
has been investigated recently. Here, we discuss the interaction between dark strings and 
standard cosmic strings and show that they can form bound states. In the presence of the 
interaction term, a Bogomolny-Prasad-Sommerfield (BPS) bound exists that depends on the 
interaction parameter and we observe that the attractive interaction between dark strings 
and cosmic strings is most efficient if the two strings are identical. Moreover, our model 
allows for dark string solution that can lower their energy by coupling to an electromagnetic 



field. We also investigate the gravitational properties of these solutions and show that they 
become supermassive with a singular space-time for values of the gravitational coupling 



larger as compared to the non-interacting case. Moreover, the deficit angle of the solutions 



decreases with increasing interaction. 

On 

in 

: 1 Introduction 

There is strong observational evidence pQ that approximately 22% of the total energy density 
of the universe is in the form of dark matter. Up until now it is unclear what this dark matter 
should be made of. One of the favourite candidates are Weakly Interacting Massive Particles 
(WIMPs) which arise in extensions of the Standard Model. Recently, new theoretical models 
of the dark matter sector have been proposed [2], in which the Standard Model is coupled 
to the dark sector via an attractive interaction term. These models have been motivated by 
new astrophysical observations [3] which show an excess in electronic production in the galaxy. 
Depending on the experiment, the energy of these excess electrons is between a few GeV and 
a few TeV. One possible explanation for these observations is the annihilation of dark matter 
into electrons. Below the GeV scale, the interaction term in these models is basically of the 
form of a direct coupling between the U(l) field strength tensor of the dark matter sector and 
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the U(l) field strength tensor of electromagnetism. The U(l) symmetry of the dark sector 
has to be spontaneously broken, otherwise a "dark photon" background leading to observable 
consequences would exist. 

Consequently, it has been shown that the dark sector can have string-like solutions, denomi- 
nated "dark strings" and the observational consequences of the interaction of these dark strings 
with the Standard Model have been discussed [I]. 

Topological defects are believed to have formed in the numerous phase transitions in the 
early universe due to the Kibble mechanism [5]. While magnetic monopoles and domain walls, 
which result from the spontaneous symmetry breaking of a spherical and parity symmetry, 
respectively, are catastrophic for the universe since they would overdose it, cosmic strings are 
an acceptable remnant from the early universe. These objects form whenever an axial symmetry 
gets spontaneously broken and (due to topological arguments) are either infinitely long or exist 
in the form of cosmic string loops. Numerical simulations of the evolution of cosmic string 
networks have shown that these reach a scaling solution, i.e. their contribution to the total 
energy density of the universe becomes constant at some stage. The main mechanism that 
allows cosmic string networks to reach this scaling solution is the formation of cosmic string 
loops due to self-intersection and the consequent decay of these loops under the emission of 
gravitational radiation. 

For some time, cosmic strings were believed to be responsible for the structure formation 
in the universe. New Cosmic Microwave background (CMB) data, however, clearly shows that 
the theoretical power spectrum associated to Cosmic strings is in stark contrast to the observed 
power spectrum. However, there has been a recent revival of cosmic strings since it is now 
believed that cosmic strings might be linked to the fundamental strings of string theory [6]. 

While perturbative fundamental strings were excluded to be observable on cosmic scales for 
many reasons [7J, there are now new theories containing extra dimensions, so-called brane world 
model, that allow to lower the fundamental Planck scale down to the TeV scale. This and the 
observation that cosmic strings generically form at the end of inflation in inflationary models 
resulting from String Theory [8] and Supersymmetric Grand Unified Theories [9] has boosted the 
interest in comic string solutions again. The interaction of cosmic strings has been investigated 
in the context of field theoretical models describing bound systems of D- and F-strings, so-called 
p-q-strings pUl E] • 

Here we study the interaction of cosmic strings with dark string solutions. On a field theoret- 
ical level, the model is similar to the one used in [10} 111), however in this paper, the interaction 
between the strings is mediated via the gauge fields (and gravity), while the strings in |10} lllj 
interact via a potential (and gravity). While the dark sector in our model is an Abelian Higgs 
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model, we should in principle couple the corresponding solutions to the electromagnetic field 
of the Standard model. Here, we assume the U(l) symmetry to be spontaneously broken and 
the corresponding Abelian-Higgs model to possess cosmic string solutions. In fact, we employ 
the strategy that the Standard model (and its semilocal brother) have string-like solutions [T2] 
which share many features with the solutions of the U(l) toy model. Solutions in this U(l) 
Abelian-Higgs model have been first discussed in [13]. These solutions have a magnetic field 
with a quantized magnetic flux and are topolgical solitons in the sense that they have a topolog- 
ical charge associated to them. When the gauge boson mass is equal to the Higgs boson mass 
in this theory, the string solutions fulfill an energy bound, the Bogomolny-Prasad-Sommerfield 
(BPS) bound [H], such that the energy per unit length is directly proportional to the topological 
charge. These solutions have been discussed extensively and their gravitational properties have 
also been investigated [15 \ fTT] . The main feature of the space-time around a cosmic string is 
that it is locally flat, but globally possesses a deficit angle 5 that is directly proportional to the 
energy per unit length of the solutions /j,: 5 = 8nGfi, where G is Newton's constant. This leads 
e.g. to gravitational lensing effects that should make it possible to detect cosmic strings in the 
universe. Interestingly, globally regular gravitating strings exist only as long as the solutions are 
not too massive. If they become too massive (or the gravitational coupling becomes too large) 
the deficit angle is larger than 2tt and the space-time is singular [18] . It has been noted [161 [T7] 
that solutions that are BPS in flat space-time are also BPS is curved space-time, i.e. fulfill the 
same energy bound. 

Our paper is organized as follows: in Section 2, we give the model, the equations and discuss 
the BPS bound. In Section 3, we present our numerical results and Section 4 contains our 
conclusions and an outlook. 



We study the interaction of a U(l) Abelian-Higgs field model, which has cosmic string solutions, 
with the low energy dark sector, which is also a U(l) Abelian-Higgs model in flat and curved 
space-time, respectively. Note that the U(l) model is a toy model here for standard model-like 
theories with gauge group SU(2)x U(l), which also contain string solutions. Examples would 
be semilocal strings in the SU(2) 9 ; b a ; x U(l) model and electroweak strings in the SU(2); oca ; x 
U(l) model [12], respectively. We believe that our toy model captures the qualitative features 
of these theories. 

The model we are studying is given by the following action: 



2 The model 




(i) 
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where R is the Ricci scalar and G denotes Newton's constant. The matter Lagrangian C m reads: 

L m = D^(D^y - l -F^ + D^(D^T ~ \h^H^ - + £ -F^H^ (2) 

with the covariant derivatives D^cf) = — ieiA^cp, = — iei^^ and the field strength 
tensors F^ v = d fM A u — d u A^, H^ v = d^a u — dya^ of the two U(l) gauge potential A^, with 
coupling constants e\ and <f> and £ are complex scalar fields (Higgs fields) with potential 

ma^^-tf+i^-^ 2 ) 2 ( 3 ) 

The term proportional to e is the interaction term [3]. To be compatible with observations, e 
should be on the order of 10~ 3 . 

In the following, we associate the dark strings to the fields A^ and (ft, while the standard 
cosmic strings are described by the fields and £. The Higgs fields have masses Mh,% = V^iVii 
while the gauge boson masses are M\y,i = V^e^, i = 1,2. 



2.1 The Ansatz 

In the following we shall analyse the system of coupled differential equations associated with the 
gravitationally coupled system described above. This system will contain the Euler-Lagrange 
equations for the matter fields and the Einstein equations for the metric fields. In order to do 
that, let us write down the matter and gravitational fields as shown below. The most general, 
cylindrically symmetric line element invariant under boosts along the z— direction is: 

ds 2 = N 2 (p)dt 2 - dp 2 - L 2 (p)d^ 2 - N 2 (p)dz 2 . (4) 

The non- vanishing components of the Ricci tensor R v ^ then read |15j : 

_ (LNN'Y 2N" L" {N 2 V)> 

°~ WIT ' "~ aT~T ' *~ WT ' Kz - Ko 

where the prime denotes the derivative with respect to p. 
For the matter and gauge fields, we have [13j : 



( f>(p, ( p)=r n h(p)e i ^ , S(p )( p) = mf(p)e mip (6) 

A.dx 11 = — (n - P(p))dip , dndx 11 = —(m - R(p))d<p . (7) 

n and m are integers indexing the vorticity of the two Higgs fields around the z— axis. 
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2.2 Equations of motion 

We define the following dimensionless variable and function: 



x = eirjtp , L[x) = rneiL(p) . (8) 

Then, the total Lagrangian C m — > C m /(rjfe 2 ) depends only on the following dimensionless 
coupling constants 

7 = 8ttG^ , A = — — = - , i = 1,2 , (9) 

and the dimensionless ratios of the coupling constants and vacuum expectation values, respec- 
tively 

9 = — , Q = — ■ (10) 
ei rji 

Varying the action with respect to the matter fields and metric functions, we obtain a system of 
six non- linear differential equations. The Euler-Lagrange equations for the matter field functions 
read: 

(N 2 LhJ P^h Pi 2 _ 
N 2 L L 2 + 2 [ ' [ ' 

= 2h 2 P + 2egRf 2 , (13) 



N 2 L 








(T) 




'N 2 R'V 


N 2 \ 





(l-^)^-^J = 2g 2 f I R + 2egPh I , (14) 

where the prime now and in the following denotes the derivative with respect to x. 
We use the Einstein equations in the following form: 

Rfiu = -7 (r^v - \a^T^j , n,v = t,x,(p,z (15) 

where T is the trace of the energy momentum tensor T = and the non-vanishing components 
of the energy-momentum tensor are (we use the notation of [15]) with i = 1, 2: 

-^o = &s ~\~ &v ~\~ &w ~l~ > T x = e s e v + e w + u 

T£ = e s -e v -e w + u , T z z = T ° (16) 



where 



(*)+(/) > ^ = nr + 2^i2---p- > e - = -^ + -L2- (17) 
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and 



We then obtain 



(LNN'Y 
N 2 L 
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(P') 2 {R 1 ? z R'P' 

~2J? + 2gH? ~ ~g~T? r 



— u 



and: 



(N 2 L 



2 T'Y 



-7 



2h 2 P 2 2R 2 f 2 {P^_ {R') 2 



eR'P' 



L 2 L 2 2L 2 2g 2 L 2 g L 2 



+ u 



N 2 L 

2.3 Boundary conditions 

The requirement of regularity at the origin leads to the following boundary conditions: 

h(0) = 0, /(0) = , P(0) = n , i?(0) = m 

for the matter fields and 

N(0) = 1, JV'(O) = 0, L(0) = , L'(0) = 1 . 
for the metric fields. The finiteness of the energy per unit length requires: 

h(oo) = 1, /(oo) = g , P(oo) = , P(oo) = . 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



2.4 Energy per unit length, magnetic fields and deficit angle 

We define as inertial energy per unit length of a solution describing the interaction of a dark 
string with winding n and a cosmic string with winding m: 



(24) 



where g% is the determinant of the 2 + 1-dimensional space-time given by (t,x,<p). This then 
reads: 

[•OO 

(25) 



poo 

M (n,m) = 27r / NL (e s + E v + E w + u) dx 

Jo 



Note that the string tension T = J \f—gz T*dxd<p is equal to the energy per unit length. In flat 
space-time (G = 0) and s = 0, the energy per unit length of the solution is given by 



(26) 



where g\ and g^ are functions that depend only weakly on j3\ and 02, respectively with <?i(2) = 1 
and 52(2) = 1 in the BPS limit. In the following, we will set 771 = 1 without loss of generality. 
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We can define the binding energy of an ra-dark string with an m-cosmic string as 

^ = fj,M _ n/i (i,o) _ mfI (o,i) (27) 



The magnetic fields associated to the solution can be given when noting that the gauge part 
of the Lagrangian density can be rewritten as follows [3J: 

- - + E -F^H^ =► —F^F^ e 2 )H^ v H^ (28) 



with = d^A u - duA^ where = - ea M . 

The magnetic fields associated to the fields A^ and have only a component in z-direction. 
These components read: 

-P'(x) + ^R'(x) , R'( x ) 

B z (x) = '.\ and b z {x) = -^^-^-^L (29 ) 

e\L{x) e2L(x) 

respectively. The corresponding magnetic fluxes J d 2 x B are 

~ 2ir ( e \ , / 7; 2-irm 

$ = — n m and p = V 1 - e 2 , (30) 

ei V g J ' e 2 

respectively. Obviously, these magnetic fluxes are not quantized for generic e. 

Finally, the deficit angle 5 = ^ttG^l of the solution can be read of directly from the derivative 

of the metric function L{x). For string-like solutions, the metric functions behave like N(x — > 

00) — > c\ and L(x — ► 00) — > c 2 x + C3, where ci, C2 and C3 are constants. The deficit angle is then 

given by: 

S = 2tt(1 - L% =00 ) = 2tt(1 - c 2 ) . (31) 
2.5 The Bogomolny-Prasad-Sommerfield (BPS) bound 

For e = 0, the model has BPS solutions which satisfy the energy bound \x = 2irn + 2itm both in 

flat [H] as well as in curved space-time [T7j for f}\ = 2, = 2. 

Here we will discuss the case f = h {j3\ = f3 2 = , P = R (n = m, g = 1) for e ^ 0. In flat 

space-time, the functions N = 1 and L = x and the corresponding BPS equations read 

P/> P' 
h' = — , (l-e)— = h 2 -l (32) 
x x 

The solutions then fulfill the energy bound E = 2irn + 27rm = 4-7rn for 

" = rb <33) 

In curved space-time, i.e. for 7 / 0, the above bound is still a BPS bound. The corresponding 
matter equations read 

h'=™ , {l-e)^ = h 2 -l. (34) 
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We have e v = u and e w = e s such that T£ = T£ = 0. Hence N(x) = 1, while L fulfills the 
following equation: 

The deficit angle then reads 

<5 = 27r 7 (n + m) . (36) 



2.6 Dark strings coupled to an unbroken U(l) symmetry 

For /?2 = 0, q = (with 771 finite), we have / = and the corresponding U(l) symmetry remains 
unbroken. For s = 0, m = 0, we would then have R = 0. However, for e ^ 0, m = the energy 
of the solution can be lowered by a non-vanishing derivative of R, i.e. R' ^ 0. This corresponds 
to a non- vanishing magnetic field in z-direction. Due to the attractive nature of the interaction, 
the dark string can thus lower its energy by coupling to a non-vanishing magnetic field inside 
the string core. The corresponding boundary conditions for R then read: 

R(Q) = and R'\ x=00 = (37) 

Note that we don't need R but only R' to vanish at infinity to fulfill the requirement of finiteness 
of the energy. (|17p then suggests that the e v part of the energy density is minimized for R' = egP' 
inside the dark string core. Using the boundary conditions (I2ip . (|37[) we find that R = eg(P-n). 
The energy of an (n, 0) string is thus smaller for e / than for e = 0. The same is, of course, 
also true for the inverse case, where n = and m ^ and R replaced by P in the boundary 
conditions ([371) . 



3 Numerical results 

In all our numerical calculations, we have set q = 1, g = 1. We have solved the set of cou- 
pled ordinary differential equations numerically using the ordinary differential equation solver 
COLSYS [20]. Numerical errors are typically on the order of 10 -8 — 10 -9 . 



3.1 7 = 

To begin with, we present our results for cosmic-dark string systems in flat space-time 7 = 0. 
In this case, we have N = 1 and L = x. 

We have first studied solutions with m = 0, = 0. As discussed above, the field / = 0, but 
R ^ 0. A typical solutions of this type is shown in FiglTJ where we present the functions P(x) 
and h(x) associated to the dark strings as well as the field R{x) associated to the unbroken U(l) 
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Figure 1: The functions of a dark string solution for re = 3, m = 0, fa = 2, fa = and e = 0.1. 
The fields associated to the dark string are P{x) and h(x), while the scalar field of the unbroken 
U(l) symmetry vanishes identically f(x) = 0. The gauge field function R(x) however is related 
to the function P{x) by R(x)/e = P{x) — re. 

symmetry. The energy per unit length of the dark strings is lowered from = 3 for e = 

to // 3 >°) = 2.99080 for e = 0.1. 

Next, we have investigated the existence of bound states in our system. Following the 
definition of the binding energy ([2*7]) . we present our results for e ^ and 0i = fa = P. In the 
e = limit, /x["; m) > for f3 > 2, < for (3 < 2 and = for (3 = 2 (the BPS limit). 

This is different here as the tables below suggest. In Tables 1,2,3 and 4,5,6 we give the binding 
energy for e = 0.001 and e = 0.01, respectively, for different choices of (3 and n and m. 



(n,m) 


1 


2 


3 


4 


5 


1 


-0.0004 


-0.0006 


-0.0008 


-0.0009 


-0.0009 


2 


-0.0006 


-0.0011 


-0.0014 


-0.0016 


-0.0017 


3 


-0.0008 


-0.0014 


-0.0028 


-0.0022 


-0.0024 


4 


-0.0009 


-0.0016 


-0.0022 


-0.0026 


-0.0030 


5 


-0.0009 


-0.0017 


-0.0024 


-0.0030 


-0.0034 



Table 1: The value of the binding energy /xL ' in units of 2n for fa = = 2, e = 0.001 and 
different choices of n and m. 
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(n,m) 


1 


2 


3 


4 


5 


1 


-0.0005 


0.0566 


0.1319 


0.2169 


0.3082 


2 


0.0566 


0.1134 


0.1886 


0.2734 


0.3646 


3 


0.1319 


0.1886 


0.2635 


0.3482 


0.4393 


4 


0.2169 


0.2734 


0.3482 


0.4328 


0.5237 


5 


0.3082 


0.3646 


0.4393 


0.5237 


0.6144 



Table 2: The value of the binding energy ir bin in units of 2ir for (3\ = Pi = 3, e = 0.001 and 
different choices of n and m. 



(n,m) 


1 


2 


3 


4 


5 


1 


-0.0005 


0.1049 


0.2448 


0.4033 


0.5738 


2 


0.1049 


0.2099 


0.3497 


0.5081 


0.6784 


3 


0.2448 


0.3497 


0.4893 


0.6475 


0.8177 


4 


0.4033 


0.5081 


0.6475 


0.8055 


0.9755 


5 


0.5738 


0.6784 


0.8177 


0.9755 


1.1454 



Table 3: The value of the binding energy fi bin in units of 2tt for (3\ = = 4, e = 0.001 and 
different choices of n and m. 

Increasing e increases the attractive interaction between the dark string and the comic string. 
For e = and P\= Pi = 2, the strings do not interact, while for the same choice of the /3s and 
£ ^ they form bound states for all values of n and m. We observe that the binding becomes 
stronger with increasing n + m. This is related to the fact that the solutions' mass increases with 
increasing winding such that the binding mechanism is more efficient. Moreover, the binding 
is stronger for an (n, n) string as compared to an (n + 1, n — 1) string and the binding is more 
effective for an (n + 1, n — 1) string than for an (n + 2, n — 2) string etc. Apparently, dark strings 
bind stronger to cosmic strings if they have the same winding. 

For e = and Pi = fa = 3 and 0\ = fii = 4, respectively, the strings repel since in this 
case the radius of the flux tube core is larger than that of the scalar core. Apparently, in this 
case the attractive nature of the new interaction (e / 0) is not strong enough to overcome the 
repulsion for these choices of (3\ = expect in the case n = m = 1. Apparently, the attractive 
interaction is just about able to overcome the repulsion that two strings would exert on each 
other for e = and > 2, i = 1, 2. 

We have next determined the critical value of Pi = fa = P = Per for which the transition 
between bound and unbound strings takes place. For e = 0, this happens at P cr {e = 0) = 2 such 
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(n,m) 


1 


2 


3 


4 


5 


1 


-0.0042 


-0.0064 


-0.007 


-0.0085 


-0.0090 


2 


-0.0064 


-0.0109 


-0.0138 


-0.0158 


-0.0171 


3 


-0.0077 


-0.0138 


-0.0183 


-0.0216 


-0.0240 


4 


-0.0085 


-0.0158 


-0.0216 


-0.0262 


-0.0297 


5 


-0.0090 


-0.0171 


-0.0240 


-0.0297 


-0.0344 



Table 4: The value of the binding energy in units of 2ir for fa = fa = 2, s = 0.01 and 

different choices of n and m. 



(n,m) 


1 


2 


3 


4 


5 


1 


-0.0046 


0.0501 


0.1240 


0.2081 


0.2988 


2 


0.0501 


0.1023 


0.1743 


0.2571 


0.3468 


3 


0.1240 


0.1743 


0.2246 


0.3258 


0.4142 


4 


0.2081 


0.2571 


0.3258 


0.4055 


0.4926 


5 


0.2988 


0.3468 


0.4142 


0.4926 


0.5784 



Table 5: The value of the binding energy fij- ' in units of 2n for fa = fa = 3, e = 0.01 and 
different choices of n and m. 

that for (5 < 2, the binding energy is negative, while for (3 > 2, the binding energy is positive. 
Here, the additional attractive interaction allows us to increase the ratio between Higgs and 
gauge boson mass to values larger than (5 = 2 before the dark string and cosmic string become 
unbound. To understand the dependence of (3 cr on e, we have chosen solutions with n + m 
constant and have determined the corresponding (3 cr for e = 0.01. Our results are given in Fig|2l 
We plot the mass ^( n,m ) of solutions with (n, m) = (4+/c, 4— £;), k = 0, 1, 2, 3. For comparison, 
we also plot (4 + k) times the mass of an (1, 0) solution plus (4 — k) times the mass of an (0, 1) 
solution. At the intersection points of this latter curve with the (4 + k, 4 — A;) curves, we have 
fihin = 0- For all solutions, we observe that the mass increases linearly with fa = fa — P in 
the range (3 S [2 : 2.1] and that the mass of an (4 + k\, 4 — k\) solution is higher than that of a 
(4 + A>2 5 4 — A>2) solution for k\ > &2- Moreover, j3 cr , i.e. the value at which the strings become 
unbound is largest for the solution with n = 4, m = 4 and decreases for k increasing. Again, 
this results from the fact that dark strings "bind best" to cosmic strings which have the same 
winding. 

The magnetic fields (f29|) change as compared to the e = case. We observe that -B^(O) 
decreases with increasing e, while b z (0) increases with increasing e. In addition the core radii of 
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(n,m) 


1 


2 


3 


4 


5 


1 


-0.0049 


0.0978 


0.2362 


0.3936 


0.5634 


2 


0.0978 


0.1978 


0.3341 


0.4900 


0.6587 


3 


0.2362 


0.3331 


0.4684 


0.6226 


0.7899 


4 


0.3936 


0.4900 


0.6226 


0.7752 


0.9410 


5 


0.5634 


0.6587 


0.7899 


0.9410 


1.1053 



Table 6: The value of the binding energy fA,' in units of 2ir for fix = Pi = 4, e = 0.01 and 
different choices of n and m. 

both flux tubes decrease with increasing e. 
3.2 7^0 

We have first studied the dependence of the deficit angle 5 on the interaction parameter e. Our 
results are shown in Fig J3] for 7 = 0. 2, /?i = /?2 = 2 and different choices of n and m. For e = 0, 
the deficit angle 5(e = 0) is given by (|36p . Here, we plot the difference between 5(e = 0) and 
S(e 7^ 0). The deficit angle decreases with increasing e and apparently decreases (approximately) 
linearly with e. The larger the sum n + m, the stronger 5 decreases with e, which is related to 
the fact that the higher n + m, the more massive the solutions are and the more effective is the 
attractive interaction between the dark string and the cosmic string. Moreover, for a fixed sum 
n + m, the deficit angle decreases stronger for an (n, n) string as compared to an (n + 1, n — 1) 
string. The reason for this is that the binding mechanism is better for a dark string and a cosmic 
string with equal winding. This has also been observed before in models describing interacting 
cosmic strings when the interaction is mediated via a potential term [10] [TT] . 

Another important feature of cosmic strings in a curved space-time is that globally regular 
gravitating solutions exist only up to a maximal value of the gravitational coupling 7 = 7 ma3; (e) . 
For 7 > jmax, the- deficit angle becomes larger than 2ir, i.e. the solutions become singular. For 
e = 0, these solutions have been denominated "supermassive" string solutions [18] . For e = 
and Pi = P2 = 2 we have 7 ma x(e = 0) = l/(n + m). For e > 0, the solutions exist on a larger 
interval of 7. We find e.g. for n = 1, m = 1 and Pi = pz = 2 that 7 max (e = 0.01) « 0.501, 
while jmax(£ = 0.1) ~ 0.515. This is not surprising since the total energy per unit length of 
the solutions decreases for increasing e. Since the deficit angle is proportional to the product 
of the gravitational coupling and the energy per unit length, the gravitational coupling can be 
increased stronger before 5 becomes equal to 2ir. 

For n / and m = 0, we find the same phenomenon as in flat space-time: dark string 
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Figure 2: The value of the inertial mass ^( 71 ' m ) in units of 2tt is given in dependence on (5\ = 
for e = 0.01 and different choices of n and m. Note that n + m = 8. For comparison, we also 
plot n times the mass of the (1, 0) solution plus m times the mass of the (0, 1) solution. At the 
intersection of the (n, m) curves with this latter curve, we have M&fe = 0. 
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Figure 3: The difference between the deficit angle S/(2tt) and the deficit angle for e = 0, 
5(e = 0)/(2tt) is given as function of e for 7 = 0.2, f3\ = fa = 2 and different choices of (n,m). 
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j2,0) Melvin 
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e 

Figure 4: The value of the inertial mass ^( n,m ) of string and Melvin solutions is shown in 
dependence on the interaction parameter e and 7 = 0.5. For n = m = 1, we have chosen 
01 = 02 = 2, while for n = 2, m = 0, we have 0\ = 2, 0% = 0. 

solutions can lower their energy by coupling to the gauge field of an unbroken U(l) symmetry. 
In FiglU we compare the masses of string-like solutions for n = 1, m = 1 and n = 2, m = 0, 
respectively, for 7 = 0.5. For n = 1, m = 1, we have 0\ = 02 = 2, while for n = 2, m = 0, 
we have /3i = 2, 02 = 0. In both cases, the mass of the solution is decreasing with increasing 
interaction. It is thus apparent that also in curved space-time, a dark string can lower its energy 
when coupling to an unbroken U(l) gauge symmetry. Moreover, the total energy of the (1, 1) 
system is always lower than that of the (2, 0) system. Hence, the binding is more efficient if a 
dark string and a cosmic string couple than if a dark string couples to a pure gauge field. 

It has been observed that string-like solutions with a behaviour of the metric functions 
N(x — > 00) — ► c\ and L(x — > 00) — ► ci% + C3 have a "shadow solution" in the form of so- 
called Melvin solutions |X5|, I16j. In contrast to string- like solutions, Melvin solutions have no 
flat space-time counterparts and their metric functions behave as N(x — * 00) — > a\x 2 ^ and 
L(x — > 00) — > a2X -1 / 3 , where a\ and a2 are constants. These solutions are also present in 
systems where cosmic strings interact via a potential term and we also find them here. In 
FigiU we give the inertial mass of these solutions for two different choices of n and m. When 
comparing to the string solutions, it is apparent that the mass of the Melvin solutions is higher 
for all choices of the interaction parameter. Moreover, in contrast to string solutions, the inertial 
mass of the Melvin solutions increases for increasing interaction. This has also been observed in 
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where two cosmic strings interact via a potential. The stronger the attractive interaction 
between the strings, the higher the inertial mass of the Melvin solutions. In addition, the mass 
increases stronger for an (1,1) system of dark strings than for an n = 2 dark string coupled 
to the gauge field of the unbroken U(l) symmetry. This is related to the observation that the 
binding is stronger for the (1, 1) system than for the (2,0) system. 

4 Conclusions and Outlook 

In this paper, we have studied the interaction of a dark string with a cosmic string where the 
interaction is mediated by a coupling term between the field strength tensors of the respective 
solutions. This type of interaction is motivated by recent models describing the dark matter 
sector. We observe that a BPS bound exists if the dark string is identical to the cosmic string. 
In fact, we find that the attractive interaction between the two strings is most efficient in this 
particular case. In addition, the attractive interaction allows for dark-cosmic strings to exist 
for larger values of the Higgs to gauge boson ratio. The deficit angle associated to the strings 
decreases for increasing interaction and globally regular string space-times exist for higher values 
of the gravitational coupling as compared to the non-interacting case. We also find dark string 
solutions that can lower their energy by coupling to a U(l) field through the attractive interaction 
term. 

The formation of bound states is of interest for the study of the evolution of string networks. 
For p-q-strings it has been observed [21] that the formation of bound states leads to an energy 
loss mechanism that is important in the evolution of string networks towards the scaling regime. 
If the dark matter sector has indeed dark strings solutions then standard cosmic strings could 
loose energy by coupling to these strings. 

It would also be interesting to understand the influence of the type of interaction studied 
here on superconducting strings [22J. In its simplest version this would be the coupling of a 
U(l)x U(l) model (describing the superconducting string with bosonic charge carriers) coupled 
to a U(l) model describing the dark string. The field strength tensor of the dark string could 
then either be coupled to the field strength tensor associated to the broken U(l) symmetry 
(similar to what has been studied in this paper) or alternatively to the field strength tensor of 
the unbroken U(l) symmetry describing the carrier field. 
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